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Measurement of the {220} lattice-plane spacing of a 8 Si crystal 
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(Dated: October 18, 2010) 

The spacing of the {220} lattice planes of a 28 Si crystal was measured by combined x-ray and 
optical interferometry to a 3.5 x 10 -9 relative accuracy. The result is cfeo = (192014712.67 ± 0.67) 
am, at 20.0 °C and Pa. This value is greater by (1.9464 ± 0.0067) x 10 _9 d 22 o than the spacing in 
natural Si, a difference which confirms quantum mechanics calculations. Subsequently, this crystal 
has been used to determine the Avogadro constant by counting the Si atoms, a key step towards a 
realization of the mass unit based on a conventional value of the Planck or the Avogadro constants. 

PACS numbers: 06.20.-f, 06.20.Jr, 61.05.cp, 07.60.Ly 



We participated in a project ended with an accurate 
measurement of the Avogadro constant, Na, by counting 
the atoms in a silicon crystal highly enriched with 28 Si 
[H, This could make it possible to realize the mass 
unit on the base of the mass of an atom Q. Since the 
molar Planck constant is known via the measurement of 
the Rydberg constant, this measurement will also lead 
to a value of the Planck constant, h, and, in turn, to a 
kilogram realization based on a conventional value of h. 

The measurement equation is Na = nM/(pa,Q) 
where M is the molar mass, p the density, Oq the cu- 
bic unit cell volume, arj the lattice parameter, and n the 
number of atoms per unit cell. The mass of the Pt-Ir 
kilogram-prototype may have changed by about 50 p,g 
since its manufacturing, in 1889 Owing to the small- 
ness of this drift, the relative uncertainty of the Na de- 
termination must not exceed 2 x 10 -8 . 

The crystal production started in 2004, with the iso- 
tope enrichment of the SiF4 gas by the Central Design 
Bureau of Machine Building in St. Petersburg; subse- 
quently, after conversion of the enriched gas into SiH4, 
a polycrystal was grown by chemical vapor deposition at 
the Institute of Chemistry of High-Purity Substances of 
the Russian Academy of Sciences in Nizhny-Novgorod. 
Eventually, a 5 kg 28 Si boule was grown and puri- 
fied by application of the float-zone technique by the 
Leibniz- Institut fur Kristallziichtung (Berlin) in 2007 Q . 
The pulling speed was so chosen as to reduce the self- 
interstitial concentration; no doping by nitrogen was ap- 
plied. The growth axis was [100]. 

The lattice parameter ao is one of the quantities re- 
quired to determine Na', the most accurate way to mea- 
sure it is by means of combined x-ray and optical inter- 
ferometry [7h10i|. Owing to the many contributions to 
the Na uncertainty, the uncertainty of this measurement 
should be reduced to 3 x 10~ 9 a . Therefore, we extended 
the measurement capabilities of an x-ray interferometer 
to many centimeters and manufactured an interferome- 
ter with unusually long analyzer. The measured value 
was also connected to the results of measurements of the 
lattice parameter of the natural Si reference-crystals of 
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FIG. 1: The combined X-ray and optical interferometer. 
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the PTB and the INRIM 

Experimental apparatus. The combined x-ray and op- 
tical interferometer is shown in Fig.Q] It consists of three 
crystal slabs, 1.20 mm thick, so cut that the {220} planes 
are orthogonal to the crystal surfaces. X-rays from a Mo 
K Q source are split by the first crystal and recombined, 
via two transmission crystals, by the third, called the 
analyzer. When the analyzer is moved along a direction 
orthogonal to the {220} planes, a periodic variation in the 
transmitted and diffracted x-ray intensities is observed, 
the period being the diffracting-plane spacing. The ana- 
lyzer movement is an extremely difficult task; it requires 
nanoradian attitude-control and picometer vibration and 
position controls. 

The analyzer embeds front and rear mirrors, so that 
its displacement and rotation can be measured by opti- 
cal interferometry; the necessary picometer and nanora- 
dian resolutions are achieved by polarization encoding 
and phase modulation. The symmetric analyzer shape 
allows it to be reversed; measurements can be carried 
out with the x rays crossing the crystal in both direc- 
tions so that surface effects, as well as others systematic 
effects, can be investigated. 

The measurement equation is ao = Vocfeo = 
\/8m\/(2n), where (I220 is the spacing of the analyzer 
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FIG. 2: Photograph of the 28 Si analyzer on its silicon support. 
The crystal temperature is measured inside the copper block. 



diffracting-planes and n is the number of x-ray fringes 
observed in a displacement of m optical fringes having 
period A/2. The laser source operates in single-mode 
and its frequency is stabilized against that of a transi- 
tion of the 127 l2 molecule. This ensures the calibration of 
the optical interferometer with a negligible uncertainty. 
To eliminate the adverse influence of the refractive index 
of air and to ensure millikclvin temperature uniformity 
and stability, the experiment is carried out in a thermo- 
vacuum chamber. 

In practice, 0I220 is determined by comparing the un- 
known period of the x-ray fringes against the known pe- 
riod of the optical ones. This is done by measuring the 
x-ray fringe phase at the start and end of a sequence of 
increasing displacements mA/2, where m = 1, 10, 100, 
1000, 3000, and 30000. To determine the fringe phase, 
the least-squares method is applied; the input data are 
about 300 samples of six fringes, with a 100 ms inte- 
gration time and a sample duration of 30 s. Each (I220 
measurement is the average of about 9 values collected in 
measurement cycles during which the analyzer is repeat- 
edly moved back and forth along the selected displace- 
ment. The visibility of the x-ray fringes approached 50% 
with a mean brilliance of 500 counts/s/mm 2 . 

To cope with the highly demanding request of accu- 
racy, the INRIM extended the crystal-displacement capa- 
bilities of its combined x-ray and optical interferometer 
to 5 cm. This magnification made more numerous effects 
visible and reproducible. In addition, it allowed wider 
crystal parts to be surveyed, thus increasing confidence 
in the crystal perfection and in the mean d.220 value. 

X-ray interferometer. In order to exploit the large 
displacement capability, the PTB manufactured an in- 
terferometer with an unusually long analyzer, which is 
shown in Fig. [2] To nullify the error due to the different 
ways the displacement is measured by the x-ray and opti- 
cal interferometers, by projecting it on the normals to the 
mirror and diffracting planes, the front and rear mirrors 
were accurately polished parallel to the {220} planes, to 
within 10 /xrad in the worst case. The analyzer move- 



FIG. 3: Survey of the diffracting-plane strain; obverse (side 
A) and reverse (side B) analyzer-configurations. The lattice 
parameter has been measured along the dashed line, which is 
the virtual extension of the laser beam. 



ment was also sensed and controlled on-line to bisect the 
angle between these normals. Eventually, the measured 
value was corrected for the remaining imperfections. 

The analyzer has been tested for the lattice perfection. 
To this end, ^220 measurements were made over 50 sub- 
sequent crystal slabs, 1 mm wide; x-ray intensities were 
recorded by means of a multianodc photomultiplicr hav- 
ing a vertical pile of eight Nal(Tl) scintillator crystals. 
Hence, the data were processed to obtain a d 2 20 value in 
each of the 50 x 8 pixels, (1 x 1.8) mm 2 each. The rela- 
tive CI220 variations are shown in Fig. [3J Fig. 2] shows the 
relevant histogram. The obverse and reverse configura- 
tions correspond to the interferometer crystals mounted 
as they were in the boule or in a reversed arrangment. 

Results. Figure [5] shows the diffracting-plane spacing 
values measured along the lines displayed in Fig.[3j which 
indicate the laser beam propagation. These values are in- 
sensitive to the Abbe error. This error is due to parasitic 
rotations combined with an offset between the trajec- 
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FIG. 4: Histogram of the observed strains; the standard de- 
viation is 5.6 nm/m. The solid line is the best fit normal 
distribution. 
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tories of the points sensed by the x rays (indicated by 
the line in Fig. [3]) and the laser beam (indicated by the 
impinging point on the front mirror). In addition to in- 
terpolating the measurement results to obtain the ^220 
values in the points having a zero Abbe offset, parasitic 
rotations were sensed and nullified on-line. The values 
in Figs. [3] and [5] are the averages of different surveys car- 
ried out between November 2009 and July 20f 0; the error 
bars indicate the standard deviations of the data. 

The top panel of Fig. [5] shows an outlier, which sur- 
vived to the averaging; it is also visible in the correspond- 
ing panel of Fig. [3] As shown in the bottom panel, it 
disappears when the analyzer is reversed. Therefore, we 
inferred that it docs not indicate an actual g?220 variation, 
but an apparent one. An explanation is a surface effect; 
the phase difference between the interfering x rays, which 
is the basic measured quantity, records the crystal sur- 
face, as well as any extraneous material on it. Though 
this phase-contrast image is weaker by orders of magni- 
tude than the lattice image, at the sensitivity level we 
are operating it could affect the measurement result. Be- 
cause of the the cupric-ion etching used to remove the 
grinding damage, the analyzer surfaces, though flat on 
the average, are quite rough; they display a texture with 
100 /im periodicity and 10 /jm peak-to-vallcy amplitude 
[l3l ] . Since we surmized that the observed outlier origi- 
nates from this texture or from a residual contaminations, 
we removed this datum from any subsequent analysis. 

Apart from this outlier, none of the measured g?220 val- 
ues was exactly re-observed in different surveys. Addi- 
tionally, no overlapping is possible between the g?220 val- 
ues obtained with the observe and reverse mounting of 
the analyzer. Among the many crystals we examined 
11 - 13[ , this is the first one displaying a c?220 profile hav- 
ing flatness and smoothness only limited by the residual 
scattering of the data. Since this scattering is larger than 
expected from the noises of the x-ray and optical interfer- 
ometers, investigations are under way to squeeze further 
the measurement resolution. 

To acquire the mean lattice parameter of the analyzer, 
the (^220 values measured along the lines displayed in 
Fig. [3] were averaged and the corrections listed in the 
Table U were taken into account. Measurements started 
in November 2009 and were repeated in May and July 
2010; the results are compared in Fig. [6] An exemplar 
error budget is given in Table HI A detailed analysis of 
corrections and of error contributions can be found in 

Hi-til, which 

give also the results of test measurements 
aimed at establishing a firm link with the (I220 values of 
natural Si crystals used as input data for the calcula- 
tion of a self-consisted set of values of physical constants 
[3] ■ Here we draw attention to a couple of points which 
deserve particular attention. 

In the first place, it must be noted that the largest cor- 
rection is due to the diffraction of the laser beam. At this 
level of accuracy, the relation A = cjv (the symbols hav- 
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FIG. 5: Diffracting-plane spacing values measured along the 
lines indicated in Fig. [3] obverse (top) and reverse (bottom) 
analyzer-configurations. 



ing the usual meanings) is not valid; energy disperses out- 
side the region in which it would be expected to remain, 
wavefronts bend, and their spacing varies from one point 
to an other and it is different from the wavelength of a 
plane wave. Fortunately, the relevant correction depends 
only on beam divergence; not on specific characteristics of 
the beam, such as the intensity profile [TBI] . In July 2010, 
a preliminary checking of this result has been made by 
replacing the fiber collimator with a new one. The beam 
divergence is the width of its angular power-spectrum 
and was measured with the aid of a converging lens. It 
changed from 0.170(8) mrad - the value relevant to all the 
previous measurements - to 0.189(9) mrad; consequently, 
the correction changed from the 7.26(65) x 10 -9 <i22o value 
given in TableQ]to 8.92(72) x 10~ 9 d 22 o- More important, 
a beam-profile change was detected, but, as the Fig. [6] 
shows, the result of the July 7-th measurement is consis- 
tent with the previous values. 

In the second place, a new systematic effect came into 
evidence. To avoid power dissipation inside the vacuum 
chamber and temperature gradients, the pointing system 
of the laser beam, including the electro-optical crystal 
for phase modulation, are outside the chamber, stiffly 
clamped to it. However, to cut off the vibrations of 
the vacuum chamber, the experiment platform inside the 
chamber rests on three O-rings. Owing to their limited 
stiffness, the hysteresis, and to the mass of the analyzer 
carriage (about 2 kg), the analyzer displacement causes 
random misalignment and a systematic tilt of about 200 
nrad/mm between the laser beam and the optical inter- 
ferometer. In this way, in addition to increasing noise 
and scattering of the measured values, the variation of 
the lengths optical paths of the laser beam through the 
interferometer optics causes a systematic error. A sepa- 
rate experiment quantified it in (1.37 ± 1.20) x 10 -9 <i220- 
At present, we either corrected the measured values (from 
09/11/04 to 10/05/05) or applied a feedforward compen- 
sation of the tilt by counter-rotating the platform (from 
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FIG. 6: Mean of the diffracting-plane spacing values measured 
along the lines indicated in Fig. [3] the dashed line is the av- 
erage. Measurements from 09/11/04 to 10/05/05 were made 
with the obverse analyzer, the following ones with the reverse 
analyzer. On 10/07/12, the divergence of the laser beam was 
changed. Up to 10/05/05 the measured values were corrected 
for the platform tilt; afterwards, it was nullified by a counter- 
rotation. 



10/05/10 onwards). Future activity is aimed at measur- 
ing and controlling on-line the platform tilt, as well as at 
reducing the sensitivity of the optical interferometer to 
jitter and tilt of the laser beam. 

Conclusions. The final value at Pa and 20.0 °C, 



(W 28 Si) = 192014712.67(67) am, 



(1) 



is the mean of the measurement results shown in Fig. 
[6l The relative uncertainty is3.5xl0~ 9 . To act with 
caution, the uncertainty of ([1]) is the mean uncertainty of 
each single measurement in Fig. [SJ not the value reduced 
by the mean. 

Point defects, mainly carbon, oxygen, and vacancies, 
strain the crystal. The TYa determination required that 
the defect-concentration differences between the samples 
used for the lattice parameter and density measurements 
is accounted for. Therefore, we registered the distance 
from the seed crystal, 299.5 mm, of the 50 mm line along 



TABLE I: Relative correction and uncertainty, in parts per 



10 , of the 10/05/05 measured CL220 value 



Contribution 


Correction 


Uncertainty 


statistics 


0.00 


0.36 


wavelength 


-0.06 


0.03 


laser beam diffraction 


7.26 


0.65 


laser beam alignment 


1.36 


0.77 


platform tilt 


-1.37 


1.20 


Abbe's error 


0.00 


1.50 


trajectory 


1.06 


0.65 


temperature 


-0.50 


2.55 


self weigh deformation 


0.81 


0.30 


aberration 


0.00 


0.50 


total 


8.56 


3.48 



which the measured c?220 values were averaged. 

After extrapolation to an impurity free crystal, the 28 Si 
lattice parameter is larger by 1.9464(67) x 10 _6 ao than 
the parameter of the natural Si crystal WASO04, simi- 
larly extrapolated to an impurity free crystal 
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The dependence on isotopic composition is a com- 
bined effect of thermodynamics and quantum mechanics 
3 13] ■ The interatomic distance minimizes the Gibb's 
free energy with respect to the cell volume. In addition to 
the elastic energy, the free energy depends on the phoiion 
energy as well as on the entropy associated with temper- 
ature. While the elastic energy sets an equilibrium dis- 
tance independent of nuclear mass, the phonon energy 
does not. Anharmonic effects imply a greater equilib- 
rium distance and they cause thermal expansion. Since 
heavier isotopes have smaller phonon energy, they set at 
a smaller distance. Entropy increases with temperature 
and has an opposite effect. At zero temperature, only 
the zero-point phonon energy survives so that 28 Si has 
the greater lattice parameter; this is a pure quantum me- 
chanical effect. When temperature increases, the lattice 
parameter difference decreases, as a consequence of the 
increasing entropy. At room temperature, the result of 
quantum mechanics calculations of the lattice parameter 
difference, 2.03 x 10 _6 a , between 28 Si and natural Si [lU 
is in excellent agreement with the value we found. 

We received funds from the European Communitys VII 
Framework Programme ERA-NET Plus grant 217257. 
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